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SUMMARY 


Improvements  in  Airborne  techniques  for  the  measurement  of 
the  earth’s  total  magnetic  field  gradient  using  optically  pumped  magnetom¬ 
eters  are  discussed.  These  techniques  involve  fast  and  accurate  measurements 
of  small  frequency  differences  anc  s~°  applicable  to  a  large  class  of  trans¬ 
ducers  which  have  for  an  output  a  frequency  varying  with  the  input  variable. 
Description  and  design  details  are  given  of  a  newly  developed  analog  in¬ 
strument  which  allows  a  continuous  read-out  of  frequency  difference  with  a 
resolution  of  better  than  0.01  Hz. 

\ 


RESUME 


L ’auteur  discute  des  ameliorations  des  techniques  aeriennes  pour 
la  mesure  du  gradient  du  champ  magnetique  total  de  la  terre  en  utilisant 
des  magnetometres  a  pompage  optique.  Ces  techniques  comportent  des 
mesures  rapides  et  precises  de  petites  differences  de  frequence  et  sont 
applicables  a  une  vaste  classe  de  transducteurs  qui  ont  comme  sortie  une 
frequence  variant  en  fonction  de  la  variable  d’entree.  L’auteur  foumit  la 
description  et  les  details  de  conception  d’un  instrument  analogique  de 
production  recente  qui  permet  une  lecture  continue  de  la  difference  de 
frequence  avec  une  precision  superieure  a  0.01  Hz. 
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THE  FREQUENCY  DIFFERENCE  TRANSLATOR 

AN  IMPROVEMENT  IN  AIRBORNE  TECHNIQUES  FOR  MEASURING 
THE  EARTH’S  MAGNETIC  FIELD  GRADIENT 


1.0  INTRODUCTION 

An  accurate  airborne  measurement  of  the  earth’s  total  magnetic  field  and  its  gradient  poses 
numerous  technical  difficulties.  To  measure  the  gradient  component  along  a  line,  instantaneous  meas¬ 
urements  are  required  at  two  proximate  locations;  then  the  difference  of  the  measurements  has  to 
be  computed  and  divided  by  the  distance  between  the  two  locations.  The  gradient  is  more  sensitive 
to  close  magnetic  sources,  while  distant  sources  affecting  both  measuring  instruments  similarly  will 
have  a  much  smaller  effect.  The  gradient,  therefore,  can  be  indicative  of  an  ore  body  near  the  earth’s 
surface  or  even  the  presence  of  a  submarine. 

Recently,  there  has  been  a  surge  of  interest  in  airborne  methods  of  measuring  the  gradient 
of  the  earth’s  total  magnetic  field.  A  considerable  amount  of  experimental  data  has  been  obtained 
from  the  magnetometer-equipped  Convair  580  and  its  predecessor  the  North  Star  aircraft,  both 
operated  by  the  Flight  Research  Laboratory  of  the  National  Aeronautical  Establishment.  Each  of 
these  aircraft  used  optically  pumped  cesium-vapour  magnetometers;  their  output  is  a  frequency 
proportional  to  the  total  magnetic  field.  The  measurement  of  the  gradient  components  is,  therefore, 
reduced  to  the  measurement  of  a  frequency  difference.  To  facilitate  the  measuring  and  processing  of 
total  magnetic  field  gradient,  an  analog  instrument  with  high  resolution  has  been  specially  developed 
by  the  Flight  Research  Laboratory. 

This  paper  discusses  a  few  presently  used  conventional  measurement  methods  and  their 
limitations  when  used  with  optically  pumped  magnetometers.  It  also  discusses  the  design  principles, 
the  advantages  of  the  newly-developed  analog  instrument  and  its  current  and  proposed  applications. 


2.0  MAGNETIC  FIELD  MEASURING  TECHNIQUES 

The  following  discussion  of  measuring  techniques  is  applicable  to  a  large  class  of  transducers 
which  have  for  an  output  a  frequency  varying  with  an  input  variable.  To  this  class  of  transducers  belong 
optically  pumped  and  nuclear  resonance  (NMR)  magnetometers  (Refs.  1,  2,  3,  4).  The  emphasis  here  is 
on  optically  pumped  cesium  vapour  magnetometers. 


It  will  be  recalled  that  the  earth’s  magnetic  field  is  a  vector  field.  The  optically  pumped 
cesium-vapour  magnetometer  measures  only  the  magnitude  of  the  earth’s  magnetic  vector  irrespective 
of  its  orientation  (some  magnetometers  incorporate  an  angular  positioning  system  aligning  the 
magnetometer  parallel  to  the  earth’s  magnetic  vector).  Thus  the  magnitude  of  the  earth’s  magnetic 
vector  forms  a  scalar  field  referred  to  as  the  earth’s  total  magnetic  field.  The  future  discussion  relates 
to  total  magnetic  fields  which  may  originate  from  various  magnetic  dipoles  which  are  summed  with  that 
of  the  earth. 

In  order  to  discuss  methods  of  measuring  the  magnetic  field  difference  it  is  necessary  to 
describe  how  the  magnetic  field  is  measured  and  the  limitations  of  the  measurement  process. 


The  output  of  a  magnetometer  is  a  frequency,  proportional  to  the  magnetic  field  0M  ,  called 
the  Larmor  frequency  fL  with 


fL  -  aK 


(1) 


where  a  ss  3.5  Hz/7*  for  cesium  and  0M  is  in  gammas  (7)  (Ref.  5). 


*  *  10'9  Tetla 
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The  time  varying  magnetic  field  0M  (t)  as  seen  by  the  airborne  magnetometer  can  be  con¬ 
sidered  as  the  sum  of  the  earth’s  field  </>,.  (t)  and  various  interfering  fields  <f>  ( t) 

J 

(&M(t)  =  +  (2) 

To  measure  the  magnetic  field  the  Larmor  frequency  has  to  be  measured.  Usually  the  Larmor  fre¬ 
quency  is  counted  over  a  certain  fixed  time  interval  with  a  digital  frequency  counter.  The  accumulated 
count  will  correspond  to  the  average  value  of  the  field  over  this  time  interval;  the  process  is  that  of 
integration  with  respect  to  time.  The  averaging  time  must  be  long  enough  to  reduce  the  effect  of  high 
frequency  interference;  but,  at  the  same  time,  it  must  be  short  enough  to  follow  the  field  changes  with 
time  such  as  those  caused  by  the  movement  of  the  aircraft  through  the  variable  earth’s  field. 

The  digital  counter  used  to  measure  the  frequency  is  one  of  the  main  building  blocks  of  the 
instrumentation  and  its  behaviour  will  be  briefly  described.  The  frequency  response  of  the  counter 
can  be  calculated  as  follows:  assume  an  instantaneous  frequency  f(t)  of  a  frequency-modulated  signal, 

f(t)  =  fQ  +  Asin27rFt  (3) 

where  A  <  f  and  F  is  the  frequency  of  modulation. 

The  average  value  fn  of  f(t)  is  corresponding  to  the  nth  sample  of  duration  T  and  it  is  defined 


■  tL 


(n+l)T  SUITtFT 

(fo  +  A  sin2?rFt)dt  =  fQ  +  A  sin[27rF(n+V6)T]  - 


It  follows  that  S(F),  the  digital  counter  response  to  a  sinusoidal  frequency  modulation  Asin2irFt  of 
the  input  frequency  is: 


S(F)  =  A 


sinjrFT 


In  order  to  estimate  the  upper  bound  of  error  contributed  by  an  interfering  magnetic  field 
assume  the  field  to  be  of  the  form: 

0j(t)  =  01sin27rF1t  (( 

From  Equation  5,  the  contribution  of  the  interference  to  the  measured  value  of  the  field  with  a 
counter  is: 

sinjrFjT  <t>{ 

0  -  <  -  (' 

1  ttFjT  ttFjT 

For  example,  to  find  the  effect  of  the  Convair  580  propeller  interference  (with  0(  =  I7 
and  F,  =  17  Hz)  one  has: 


ttF,T  53T 

Thus,  for  a  O.Oly  maximum  contribution  to  0M  (t),  T  would  have  to  be  approximately  2  seconds; 
similarly,  for  a  O.I7  contribution,  a  T  of  0.2  seconds  would  be  sufficient.  It  follows  that,  if  the  magni¬ 
tude  and  the  frequency  of  an  interference  are  known,  an  averaging  period  T  can  be  selected  to  reduce 
its  contribution  below  a  specified  level. 
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The  uncertainty  of  count  inherent  in  an  electronic  digital  counter  is  one  count.  With  a  gate 
or  measuring  period  T  and  average  input  frequency  f  the  counter  will  accumulate  N  counts  according 
to: 

N  =  f  T  (9) 

For  a  one  count  uncertainty,  the  frequency  uncertainty  is 

Af  =  ~  (10) 

T 

If  the  input  frequency  is  multiplied  by  a  factor  K  then  the  frequency  uncertainty  relative  to  the  input 
becomes 


Afm  KT 


The  corresponding  magnetic  field  resolution  is 

Af.  i 

in  X 

A0 - - -  (12) 

a  aKT 

where  a  is  defined  by  Equation  (1). 

Thus,  an  increase  in  the  resolution  of  frequency  measurements  may  be  obtained  by  multi¬ 
plying  the  input  frequency  by  a  factor  K.  For  example,  to  obtain  O.I7  resolution  with  a  0.1  second 
gate  time  T,  a  K  of  approximately  29  is  required. 

The  frequency  multiplication  is  achieved,  in  practice,  with  a  phase  locked  loop,  referred  to 
later  as  PLL.  There  is  a  practical  upper  limit  to  the  multiplication  factor  K  that  can  be  used.  This 
limit  depends  on  the  phase  noise  and  disturbances  added  by  the  PLL  and  on  the  phase  noise  of  the 
input  frequency  (Ref.  8). 

Small  frequency  changes  can  be  measured  incrementally  with  a  high  resolution  frequency- 
to-voltage  converter  (Ref.  6)  such  as  a  delay  line  frequency  discriminator,  later  referred  to  as  DLD. 

In  a  DLD  a  periodic  input,  such  as  V(t)  =  Asin27rft  is  delayed  by  r  by  means  of  a  delay  line  to  form 
V(t-r)  =  Asin27rf(t-T).  The  delayed  input  and  the  non-delayed  input  are  connected  to  a  phase 
detector.  The  output  of  a  phase  detector  is  a  voltage  vq  proportional  to  the  phase  difference  tp  be¬ 
tween  the  two  inputs.  The  phase  difference  is  given  by: 

\p  =  avo  =  27rft-  27rf(t-r)  =  2irfr  (13) 

The  phase  is  linear  with  frequency  and  is  a  periodic  function  of  angle  with  a  period  of  2ir, 
therefore,  to  each  2ir  phase  change  there  corresponds  a  frequency  change  Af  such  that  Af  •  r  =  1  or 
1 

Af  =  —  .  Equation  (13)  becomes 
r 

f 

\p  =  2?rfT  =  2tt  —  =  2w(n +0)  =  2ir6  (14) 

where  n  is  an  integer  and  6  is  a  fraction  of  Af. 

A  diagram  of  the  DLD  is  shown  in  Figure  1  A.  The  input  frequency  is  converted  by  a  pulse  shaper 
into  a  train  of  short  pulses  corresponding  to,  say,  positive  zero  crossings.  These  pulses  are  then  delayed 
in  a  delay  line.  The  delayed  and  non-delayed  pulses  are  used  to  set  and  reset  a  flip-flop  operating  as  a 


> 
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conventional  phase  detector.  Its  duty  ratio  is  directly  proportional  to  the  phase  difference  between 
the  two  trains  of  pulses.  The  output  of  the  flip-flop  operates  a  precision  switch,  which  switches 
between  two  reference  levels:  ground  and  VR .  The  output  of  the  switch  is  filtered  with  a  low-pass 
filter.  Such  a  filter  will  have  an  effective  time  constant  much  larger  than  the  period  of  the  input 
frequency  and  it  will  provide  an  output  vQ  proportional  to  the  time  average  of  the  input.  The  output 
voltage  is  described  by  Equation  (15)  and  is  shown  in  Figure  IB: 


=  V„ 


t  +t  ,, 

on  off 


"  VR* 


(15) 


The  resolution  of  the  DLD  is  very  high,  limited  only  by  the  circuit  noise.  For  example, 
with  a  delay  of  250  microseconds  the  Af  is  4000  Hz.  By  multiplying  the  input  frequency  by  a  factor 
of  4,  a  Af  equivalent  to  1000  Hz  is  obtained.  A  20  volt  peak-to-peak  output  from  the  phase  detector 
will  then  give  a  scale  factor  of  10  mV/Hz,  permitting  the  display  of  1  Hz  as  full  scale  on  a  strip  chart 
recorder. 


The  linearity  and  repeatability  of  a  prototype  system  designed  at  the  Flight  Research  Labora¬ 
tory  was  found  to  be  better  than  0.1%.  The  design  details  of  such  a  delay  line  frequency  discriminator 
(DLD)  are  given  in  Appendix  A. 


3.0  TOTAL  MAGNETIC  FIELD  GRADIENT  MEASURING  TECHNIQUES 

The  magnitude  of  earth’s  total  magnetic  field  is  a  scalar  field;  its  gradient  is  a  vector.  To 
determine  the  gradient  uniquely  it  is  necessary  to  know  its  position,  its  magnitude  and  its  direction. 
The  co-ordinates  of  the  gradient  vector’s  origin  are  longitude,  latitude  and  altitude  obtained  from 
navigational  or  position  finding  data.  The  gradient  can  be  considered  to  be  a  vectorial  sum  of  its 
components.  Each  component,  which  is  a  spatial  partial  derivative,  can  be  approximated  by  a  simul¬ 
taneous  differential  measurement  of  the  magnetic  field  at  two  proximate  positions  located  on  a  line 
in  a  desired  direction.  The  minimum  separation  of  the  measuring  sensors  is  dictated  by  the  accuracy 
desired  for  the  gradient  and  by  the  resolution  with  which  a  magnetic  field  difference  can  be  measured. 
For  example,  on  an  aircraft  the  sensors  can  be  mounted  on  wing  tips  or  on  the  top  and  on  the  bottom 
of  the  vertical  tail  fin.  Often,  it  is  only  required  to  measure  one  component  e.g.  the  vertical,  common¬ 
ly  referred  to  as  the  “vertical  gradient”.  The  following  discussion  will  treat  the  measurement  of  only 
one  component  of  the  gradient,  with  the  understanding  that  the  same  measuring  techniques  can  be 
applied  to  other  components. 

As  explained  earlier,  the  measurement  of  the  magnetic  field  difference  with  optically 
pumped  magnetometers  is  reduced  to  the  measurement  of  frequency  differences. 

It  should  be  mentioned  that,  usually,  only  small  frequency  differences  are  measured  be¬ 
cause  the  magnetometers  are  subjected  to  similar  but  slightly  unequal  inputs. 

The  principle  of  a  typical  digital  frequency  difference  measuring  circuit  is  shown  in 
Figure  2.  In  order  to  reduce  operating  frequencies  for  the  PLL’s  and  for  the  counters,  a  down- 
conversion  can  be  used  as  shown  in  Figure  3.  Both  channels  for  fL  and  f2  are  subject  to  identical 
processing. 

Referring  to  Figure  3,  the  output  Dn  from  the  subtracting  circuit  will  be: 

_  /*  (n+l)T  /•(n+DT  /*(n+l)T 

D„  D[fj  (t)  -  fjdt  -  jfT  K[f2(t)  -  fjdt  =  KjfT  (f,(t)-f2(t)]dt 


=  KT(f,  -  f2) 


(16) 


The  gating  period  T  and  multiplication  factor  K  will  be  selected  to  satisfy  a  given  inter¬ 
ference  rejection  and  a  resolution  condition.  At  the  same  time  T  has  also  to  be  chosen  to  achieve  a 
desirable  sampling  rate. 

Reference  9  describes  an  airborne  earth’s  total  field  measuring  system  which  evolved  into 
an  airborne  gradiometer  described  in  Reference  10.  The  gradiometer  design  is  based  on  the  frequency 
difference  measuring  system,  as  shown  in  Figure  3. 

The  accuracy  to  which  a  magnetic  field  can  be  measured  will  depend  on  interference  level, 
system  noise,  timing  accuracy  and  the  ambiguity  of  one  count  in  the  digital  counter. 

An  analog  frequency  difference  measuring  system  is  shown  in  Figure  4.  It  consists  of  two 
DLD’s,  the  outputs  of  which  are  subtracted  by  means  of  two  operational  amplifiers.  A  scale  factor 
correction  is  required  to  compensate  for  different  delays  in  each  delay  line.  The  accuracy  of  such  a 
circuit  depends,  to  a  large  extent,  on  the  quality  of  the  DLD’s  used. 

Another  frequency  difference  measuring  system  is  shown  in  Figure  5.  This  arrangement 
consists  of  a  mixer  circuit  with  a  low-pass  filter  to  provide  the  frequency  difference,  followed  by  a 
digital  measuring  system  in  the  form  of  a  universal  counter-timer  (e.g.  Philips  PM6650  or  similar) 
capable  of  measuring  frequency  and  period  or  multiple  periods.  According  to  the  identity 

cos(cjjt  +  )cos  <o2t  =  VfclcosIXtOj  +  co2)t  +  *P\  1  +  cosftcjj  -  co2)t  +  l//j  ] }  (17) 

the  output  of  the  mixer  will  contain  the  sum  and  the  difference  of  cjj  and  cj2  .  After  a  low-pass  filter 
the  mixer  output  VQ  will  be: 

VQ  =  lA  cos[(tL>j  -  a>2)t  +  ]  =  Vi  cos(27rAft  +  vi'j )  (18) 

It  can  be  observed  that  the  mixing  process  preserves  the  phase  \p  j ,  which  may  have  a  large  random  com¬ 
ponent  i//(t)  (e.g.  with  a  Gaussian  distribution).  The  frequency  difference,  or  beat  frequency,  is  usually 
small,  and  to  measure  it  directly  may  require  too  much  time.  For  example,  to  measure  directly  1  Hz  to 
a  0.1  Hz  accuracy,  it  would  take  at  least  10  seconds,  severely  restricting  the  sampling  rates  available. 
Furthermore,  with  this  arrangement  multiplication  by  a  factor  K  to  increase  the  resolution  in  a  given 

Af  max 

measuring  period  T  is  difficult  when  the  frequency  ratio - is  large.  Also  a  large  ratio  complicates 

Af  min 

the  design  of  the  PLL.  Another  difficulty  with  this  arrangement  arises  from  the  ambiguity  which 
exists  as  to  whether  fj  >  f2  or  f2  >  f j ;  thus,  the  Af  is  “unsigned”. 

Period  measurements  present  similar  difficulties.  It  takes  a  long  time  to  measure  periods  of 
low  frequencies  because  the  period  tends  to  infinity  as  the  frequency  approaches  zero.  The  phase 
noise  \jj(t)  makes  the  period  readings  variable.  To  reduce  this  variability  an  average  of  N  periods  must 
be  taken,  further  increasing  the  time  of  measurement  and  adversely  affecting  the  sampling  rate. 

To  overcome  these  difficulties,  it  is  desirable  to  offset  Af  by  a  reference  frequency  fR , 
achieved,  for  example  by  adding  fR  to  f}  and  forming  (fj  +  fR  )  -  f2  =  fR  +  Af. 

The  frequency  offset  or  up-conversion  of  Af  by  fR  is  ultimately  achieved  by  the  multi¬ 
plication  (mixing)  of  two  frequencies,  a  process  which  does  preserve  the  phase  angle  \[/  of  Af  (see 
Eq.  (17)),  provided  fR  has  a  negligible  phase  noise.  Assuming  Af  «  fR  ,  during  a  measuring  period  T, 
there  are  N  phase  samples,  where  N  ss  TfR  .  Selecting  fR  and  T,  as  large  as  practical,  the  number  of 
phase  samples  N  can  be  optimized  to  obtain  phase  averages  of  a  large  number  of  samples  in  a  reason¬ 
able  measuring  period  T.  A  specialized  circuit,  referred  to  as  a  frequency  difference  translator  (FDT) 
has  been  developed  by  the  author  at  the  Flight  Research  Laboratory  for  this  purpose.  The  fR  was 
selected  to  be  111.111  kHz  although  it  could  have  been  selected  anywhere  between  100  kHz  to 
200  kHz.  Choosing  such  a  high  frequency  offset  makes  the  output  compatible  with  existing  instru¬ 
mentation  designed  to  operate  in  the  100  kHz  to  200  kHz  range,  (corresponding  to  the  cesium- 
vapour  magnetometer  output  in  the  earth’s  magnetic  field). 


4.0  FREQUENCY  DIFFERENCE  TRANSLATOR 


The  block  diagram  of  the  frequency  difference  translator  is  shown  in  Figure  6  and  the 
design  details  are  given  in  Appendix  B. 

The  principle  of  operation  is  that  of  subtracting  the  output  frequencies  of  two  phase  lock 
translation  loops.  The  output  frequency  f  can  be  expressed  as  follows: 


The  phase  lock  translation  loops  used  can  be  easily  modified  to  generate  mfj  +  fR  as  shown  in  Fig¬ 
ure  7,  resulting  in  the  frequency  difference  translator  having  an  output: 

f  =  Bfp  +  mAf  (20) 

O  tv 

A  measuring  system  is  made  up  of  a  frequency  difference  translator  supplied  with  a  stable 
reference  frequency  fR  followed  by  a  delay  line  frequency  discriminator  and  a  low-pass  filter,  the 
output  of  which  can  be  displayed  on  a  strip  chart  recorder  or  digitized  for  further  processing.  Such 
a  measuring  system  is  shown  in  Figure  8. 

To  measure  a  frequency  difference  L  -  f  j ,  exceeding  the  linear  frequency  span  or  period 
of  the  DLD’s  output,a  frequency  counter  may  be  added  (the  DLD’s  output  period  is  shown  in  the 
Figure  IB  and  denoted  Af,  not  to  be  confused  with  Af  =  f,  -  f,  used  presently).  Such  a  counter  will 
have  to  measure  the  approximate  output  frequency,  f  to  determine  by  how  many  DLD  frequency 
periods  it  differs  from  kfR  (Eq.  (19)). 

The  accuracy  of  the  system  is  limited  only  by  the  noise  generated  in  associated  mixers  and 
oscillators.  In  the  prototype  design  a  resolution  of  the  order  of  0.01  Hz  in  a  5  Hz  bandwidth  has  been 
achieved  and  a  further  improvement  in  resolution  to  0.001  Hz  in  reduced  bandwidth  can  be  antici¬ 
pated. 


The  advantages  of  such  a  system  are  mainly  stability,  high  sensitivity,  simplicity  and  speed 
of  measurement  not  previously  attainable.  It  is  essentially  a  system  operating  in  real-time  in  contrast 
to  the  system  using  digital  counters  which  provides  the  data  with  a  time  delay. 

One  of  the  important  differences  between  the  FDT  followed  by  a  DLD,  and  the  system 
operating  on  the  principle  of  digital  counting  and  subtracting,  is  that  the  former  system  is  linear  in 
frequency  and  continuous  while  the  latter  is  a  sampled  data  system.  The  output  of  the  DLD  is  a 
voltage  proportional  to  the  input  frequency.  This  voltage  can  be  processed  with  a  filter  to  suit  partic¬ 
ular  requirements.  For  example,  a  high  order  low-pass  filter  will  provide  better  rejection  of  known 
high  frequency  interferences  than  a  digital  counter  (see  Eq.  (5)).  Such  a  filter  will  reduce  the  un¬ 
wanted  effects  of  frequency  modulation  of  the  input  frequency  present  in  the  DLD’s  output. 


5.0  CURRENT  AND  PROPOSED  APPLICATIONS  OF  THE  FREQUENCY  DIFFERENCE 
TRANSLATOR 

The  most  obvious  application  is  that  of  measuring  the  earth’s  toted  magnetic  field  gradient 
and  of  variables  associated  with  the  gradient.  This  instrumentation  was  developed  primarily  for  this 
application. 

There  is  another  class  of  FDT  applications,  not  immediately  obvious,  resulting  from  the 
elimination  (subtraction)  of  magnetic  disturbances  originating  from  far  away  sources  such  as  magnetic 
micropulsations,  which  are  time  varying,  quasi-random  changes  in  the  earth’s  magnetic  field  and  which 
can  be  of  several  gammas  in  magnitude.  This  is  of  paramount  importance  in  obtaining  a  magnetically 


quiet  environment  for  experimental  purposes.  In  such  applications,  one  magnetometer  is  subjected 
to  an  external  stimulus  while  the  second  one,  located  as  close  as  possible  to  the  first  one,  acts  as  a 
“base”  or  reference.  Both  magnetometers  are  subject  to  the  same  external  “noise”  which  is  eliminated 
by  subtraction,  allowing  for  better  detectability  of  small  magnetic  changes.  With  a  single  magnetom¬ 
eter,  a  response  to  a  small  change  in  the  magnetic  field  can  be  buried  in  the  “noise”.  It  should  be 
mentioned  that  sensitive  magnetic  measurements  of  this  kind  should  be  made  in  a  nonmagnetic 
environment  removed  from  man-made  magnetic  interferences  (which  may  be  from  a  close-by  source, 
thus  generating  large  gradients). 

Magnetometer  stability  can  be  measured  by  placing  two  similar  magnetometers  a  few  feet 
apart  and  recording  the  time  history  of  the  frequency  differences.  An  assumption  is  made  that  their 
variations  are  random  and  uncorrelated  and  that  both  contribute  equally  to  the  measurement. 

Temperature  sensitivity  can  be  measured  by  subjecting  one  of  the  two  magnetometers  to 
temperature  variations,  while  the  other  is  kept  at  constant  temperature. 

Turn-on  retrace  can  be  studied  by  turning  one  magnetometer  on  and  off  periodically  and 
observing  whether  or  not  the  output  returns  to  the  same  value. 

Also,  variations  of  the  output  as  a  function  of  the  magnetometer  position  can  be  measured 
with  respect  to  the  earth’s  magnetic  field  vector.  Figure  9  shows  measured  changes  of  a  magnetom¬ 
eter  output  as  a  function  of  rotation  about  its  own  axis.  Figure  10  shows  variations  in  the  output  as 
the  magnetometer  is  tumbled  in  a  vertical  plane  parallel  to  the  North-South  bearing.  Two  “working” 
sectors  are  clearly  visible,  that  is,  the  angular  sectors  to  which  correspond  small  output  changes. 

Figure  11  shows  the  variation  when  the  magnetometer  is  rotated  in  a  horizontal  plane.  Note  that  the 
data  were  obtained  from  a  compensated,  “very  good”  magnetometer.  Other  magnetometers  may 
exhibit  much  larger  variations.  The  above  mentioned  measurements  will  allow  assessment  of  the 
possible  system  accuracy  using  a  given  magnetometer.  Furthermore,  some  corrections  can  be  intro¬ 
duced  once  a  model  is  known.  This  can  be  particularly  useful  in  a  strapped-down  configuration; 
that  is,  with  the  magnetometer  in  a  fixed  position  with  respect  to  the  aircraft,  rather  than  with  the  self 
orienting  type,  (Ref.  11). 

The  time  required  to  make  the  measurements  in  Figures  9,  10  and  11  is  relatively  short,  of 
the  order  of  a  few  minutes;  furthermore,  the  measurements  repeated  several  times  showed  a  very 
good  consistency  of  the  results.  Previous  attempts  at  similar  measurements,  using  a  counter-timer 
connected  as  shown  in  Figure  5,  resulted  in  lengthy  procedures  requiring  considerable  averaging  time 
for  each  data  point  and  the  repeated  measurements  showed  considerable  scatter. 

The  frequency  difference  translator  could  be  used  for  the  magnetic  clean-up  of  the  mag¬ 
netometer  carrying  vehicle,  in  our  case,  an  aircraft.  Assorted  electrical  systems  and  magnetic  moving 
inclusions  generate  a  man-made  magnetic  noise  which  can  reduce  the  accuracy  of  measurements; 
therefore,  it  is  essential  to  reduce  this  noise  to  a  minimum. 

The  clean-up  procedure  consists  of  placing  one  magnetometer  away  from  the  aircraft 
while  the  other  is  placed  at  various  locations  in  or  on  the  aircraft.  The  aircraft  power  is  then  turned  on 
and  various  electrical  systems  are  energized.  The  interference  from  these  systems  can  now  be  ob¬ 
served  and  minimized  by  twisting  power  leads,  using  a  “one  only”  grounding  point,  avoiding  current 
loops,  using  magnetic  shielding  etc.  Similarly,  magnetic  inclusions  in  moving  parts,  e.g.  magnetic 
steel  bolts  in  ailerons,  can  be  detected  and  when  possible  replaced  by  non-magnetic  equivalents. 

In  general,  the  knowledge  of  an  objectionable  interference  source  allows  for  a  suitable 
counteraction.  The  aim  is  to  have  in  the  aircraft  only  fixed  and  constant  magnetic  moments  whose 
behaviour  can  be  modelled  mathematically.  Use  of  the  frequency  difference  translator  allows  meas¬ 
urement  of  influences,  which  may  later  be  eliminated,  leading  ultimately  to  a  magnetically  quiet 
aircraft. 


6.0  CONCLUSIONS 


P 


The  techniques  used  to  measure  total  magnetic  field  gradient  using  optically  pumped  cesium 
vapour  magnetometers  have  been  discussed.  The  discussion  has  focussed  on  the  processing  of  the  out¬ 
put  frequency  differences.  There  are  other  transducers  and  magnetometers  with  an  output  of  fre¬ 
quency  proportional  to  the  input  variable,  to  which  the  same  processing  techniques  will  apply. 

There  are  two  basic  methods  of  measuring  frequency:  one  relying  on  counting,  the  other 
on  frequency-to-voltage  conversion.  With  the  counting  method  the  output  is  time-delayed,  since  it 
always  takes  time  to  count  and  is  only  available  in  sampled  form,  while  with  the  frequency-to- 
voltage  conversion  method  the  output  is  continuous  and  allows  for  real-time  processing. 

To  facilitate  the  measurement  of  small  frequency  differences  the  Flight  Research  Labora¬ 
tory  has  developed  a  frequency  difference  translator  which,  with  the  delay  line  frequency  discrimi¬ 
nator,  makes  a  very  versatile  instrument  providing  precise  readings.  The  output  of  such  a  system  is  a 
voltage  which  is  linear  with  respect  to  frequency  or  frequency  difference.  This  voltage  can  be  pro¬ 
cessed  with  a  frequency  filter.  The  filter  can  be  active,  of  high  order,  and  can  be  designed  to  accom¬ 
modate  specific  requirements  such  as  rejection  of  known  interferences. 


The  wide  selection  of  readily  available  “standard”  designs  for  active  frequency  filters  and 
the  precision  in  forming  a  frequency  difference  are  the  two  main  advantages  of  the  frequency  dif¬ 
ference  translator  system  over  all  the  other  systems  discussed. 

The  frequency  difference  translator  should  find  many  applications  e.g.  in  frequency  com¬ 
parison,  in  frequency  deviation  monitoring  and  in  providing  continuous  differential  readings  of  other 
transducers  with  a  variable  frequency  for  an  output.  At  the  time  of  this  writing,  the  frequency  dif¬ 
ference  translator  was  a  new  development  and  many  of  these  applications  have  yet  to  be  explored. 
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OUTPUT  OF  DELAY  LINE  DISCRIMINATOR 


COUNTER 


FIG.  3:  DIGITAL  FREQUENCY  DIFFERENCE  MEASURING  CIRCUIT  WITH  DOWN  CONVERSION 


FIG.  4:  ANALOG  FREQUENCY  DIFFERENCE  MEASURING  CIRCUIT 


TRANSLATION  LOOP  #  I 


FIG.  6:  FREQUENCY  DIFFERENCE  TRANSLATOR 


ROTATION  IN  5°  INCREMENTS  (RELATIVE  TO  45°) 


MAGNETOMETER  OUTPUT  AS  A  FUNCTION  OF  ROTATION  IN  VERTICAL  PLANE 


APPENDIX  A 


DELAY  LINE  FREQUENCY  DISCRIMINATOR  DESIGN  CONSIDERATIONS 


The  principle  of  operation  of  the  DLD  has  been  described  previously  in  the  text;  this  note 
serves  to  give  some  engineering  information  which  would  be  useful  in  designing  a  similar  device. 

A  simplified  schematic  of  a  DLD  developed  by  the  Flight  Research  Laboratory  is  shown  in 
Figure  Al.  The  phase  lock  loop  (PLL)  at  the  input  is  constructed  from  a  phase  frequency  detector 
(MC4344),  a  dual  voltage-controlled  multivibrator  (MC4324)  and  a  divided-by-four  circuit 
(SN54LS73).  Technical  data  for  similar  PLL  designs  is  contained  in  Reference  A-l. 

The  “one-shot”  (O.S.)  or  monostable  multivibrator  (SN5412I)  is  followed  by  a  buffer  B  to 
drive  the  250  ps  delay  line.  The  buffer  consists  of  four  sections  of  SN54128  wired  in  parallel.  The 
pulse  width  from  the  “one-shot”  (100  nS)  and  buffer  will  depend  on  the  bandwidth  of  the  delay  line 
used.  The  delay  line  used  was  an  ultrasonic  type  made  of  glass,  with  a  very  low  temperature  coeffi¬ 
cient  of  delay.  The  output  of  the  delay  line  is  processed  with  a  high  speed  comparator  C  (LM161). 

The  combination  of  inverters  and  NAND  gates  forms  an  edge-triggered  flip-flop  part  of  a 
conventional  phase  detector.  The  flip-flop  is  set  by  the  delayed  pulse  and  reset  by  the  non-delayed 
pulse.  The  switch  SW  and  the  Exclusive-OR  gate  in  the  reset  path  provide  for  flip-flop  reset  with 
either  the  leading  edge  or  the  trailing  edge  of  the  input  square-wave,  altering  the  flip-flop’s  duty  ratio 
by  one  half  of  a  period.  The  outputs  of  the  flip-flop  after  level  translation  circuits  L  (SN  5407), 
operate  precision  switches  S  (CD  4050),  switching  between  ground  and  a  reference  voltage  level  VR 
The  operational  amplifier  A,  with  its  associated  low-pass  filters  subtract  the  average  values  of  the 
complementary  wavetrains  generated  by  the  flip-flop  and  switches  S.  The  output  of  A,  is  a  voltage 
proportional  to  the  input  frequency  (phase),  repeating  every  Af,  which  in  this  design  is  1000  Hz. 

Note  that  the  switch  SW  changes  the  output  voltage  by  an  amount  equivalent  to  500  Hz. 

The  operational  amplifier  A2  is  also  used  as  a  low-pass  filter.  While  A!  has  a  cut-off  fre¬ 
quency  of  the  order  of  a  few  kHz,  the  response  of  A2  is  usually  tailored  to  specific  requirements, 
e.g.  60  Hz  interference  rejection. 

The  potentiometer  P  is  useful  to  off-set  the  output;  for  example,  it  may  be  required  to  set 
the  output  voltage  close  to  zero  in  order  to  display  the  output  at  a  high  sensitivity  on  a  strip-chart 
recorder. 


The  glass  delay  lines  are  available  commercially  with  delays  from  a  few  microseconds  to 
4  milliseconds.  These  delay  lines,  when  used  in  a  delay  line  frequency  discriminator,  will  provide  a 
wide  range  of  Af’s,  which  can  be  further  extended  by  either  multiplying  or  dividing  the  input 
frequency. 

The  limiting  factor  is  the  speed  with  which  the  switch  S  will  operate.  In  the  prototype 
design,  the  maximum  frequency  of  operation  is  of  the  order  of  1  MHz. 
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FIG.  A1:  SIMPLIFIED  DIAGRAM  OF  DELAY  LINE  FREQUENCY  DISCRIMINATOR 


APPENDIX  B 


FREQUENCY  DIFFERENCE  TRANSLATOR 


A  frequency  difference  translator  is  built  from  two  basic  translation  loops  as  shown  in 
Figure  6  of  the  main  text.  A  simplified  schematic  of  a  translation  loop  is  shown  in  Figure  Bl. 


There  are  two  critical  components  in  the  translation  loop  on  which  the  success  of  operation 
depends.  These  are  the  low-noise  voltage  controlled  oscillator  or  VCO,  and  a  mixer.  Careful  selection 
of  mixer  frequencies  is  also  important.  A  reasonably  high  ratio  of  fo  to  fo  -  f2  has  been  selected  to 

reduce  higher  frequency  products  coinciding  with  the  desired  output  fQ  -  f2  (Ref.  B-l).  With  the 

selected  frequencies  the  interfering  products  are  of  9th  order  or  higher.  The  selected  frequencies  are 
100  kHz  <  fj  <  200  kHz,  f2  =  1  MHz  and  1.1  MHz  <  fQ  <  1.2  MHz  for  loop  1;  and  f2  =  1.111  MHz 

and  1.211  MHz  <  f  <  1.311  MHz  for  loop  2.  One  of  the  essential  features  of  the  oscillator  is  its  low 
phase  noise.  Cj ,  C2 ,  Dt ,  D2  and  R(  form  a  voltage  doubler  in  an  auxiliary  feedback  loop  which 

biases  the  FET  transistor  into  a  linear  operating  range.  The  low-pass  filter  which  follows  the  mixer 
(MC  1596)  is  a  fifth  order  Butterworth  with  an  approximate  cut-off  freqency  of  200  kHz. 

The  output  of  the  translation  loop  is  fQ  =  fj  +  f2 .  The  f2 ’s  for  the  two  translation  loops 
10  MHz  10  MHz 

have  been  selected  to  be  - and  - respectively.  For  the  10  MHz  source  any  tempera- 

10  9 

ture  compensated  crystal  oscillator  (TCXO)  can  be  used.  Its  desired  short  time  frequency  stability 
approximately  1  X  10'7 .  The  division  is  achieved  by  conventional  digital  transistor-transistor 
(TTL).  The  subtraction  circuit,  between  the  two  translation  loop  outputs,  as  shown  in  the  Ft/.-i'  *  6 
of  the  main  text  is  achieved  with  a  mixer  (MC  1596  or  LM  1596)  followed  by  a  fifth  ord*-r  jow-pass 
Butterworth  filter,  similar  to  that  used  in  the  translation  loop. 
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G.  B1:  SIMPLIFIED  DIAGRAM  OF  FREQUENCY  TRANSLATION  LOOP 


